Structural transitions and transport-half-metallic ferromagnetism in LaMnOs at 

elevated pressure 
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By means of hybrid density functional theory we investigate the evolution of the structural, elec- 
tronic and magnetic properties of the colossal magnetoresistance (CMR) parent compound LaMn03 
under pressure. We predict a transition from a low pressure antiferromagnetic (AFM) insulator to 
a high pressure ferromagnetic (FM) transport half-metal (tHM), characterized by a large spin po- 
larization (f» 80-90 %). The FM-tHM transition is associated with a progressive quenching of the 
cooperative Jahn- Teller (JT) distortions which transform the Pnma orthorhombic phase into a per- 
fect cubic one (through a mixed phase in which JT-distorted and regular Mn06 octahedra coexist), 
and with a high-spin (S — 2, mjv/n=3.7 /^b) to low-spin (S = 1, mMn=1.7 ^b) magnetic moment col- 
lapse. These results interpret the progression of the experimentally observed non-Mott metalization 
process and open up the possibility of realizing CMR behaviors in a stoichiometric manganite. 

PACS numbers: 71.30.+h, 62.50.-p, 75.47.Gk, 72.25.-b 



I. INTRODUCTION 

Half-metallic ferromagnets (HMFs) are magnetic com- 
pounds that are simultaneously metallic and insulator, 
depending on the spin channel^ 2 - Their ability to pro- 
vide fully spin-polarized currents make them ideal candi- 
dates for spintronic applications^^ A prominent class of 
HMFs is represented by strongly correlated CMR man- 
ganites (Lai-xA^MnOs, A~Ca, Sr, or Ba)£ hole-doping 
LaMnO,3 through the substitution of La with A=C&, Sr, 
or Ba creates itinerant holes that progressively lead to 
an antiferromagnetic (AFM)-insulator to ferromagnetic 
(FM)-metal transition and critically determine the coex- 
istence of half-metallic spin imbalance and the so called 
CMR effect, i.e. a dramatic change of the electrical re- 
sistance in the presence of a magnetic field<&£ Soon af- 
ter the pioneering observation of CMR effect by Jonkcr 
and van Santen^ Zener— proposed the double exchange 
(DEX) mechanism to explain the AFM- insulator [x = 0) 
to FM-metal (0.2 < x <0.5) transition, relaying on the 
2 ~-mediatcd transfer of an electron between inequiva- 
lent Mn 3+ and Mn 4+ sites. The DEX mechanism quali- 
tatively explained the ferromagnetic interactions and the 
observed metallic behavior below the Curie temperature 
Tc, but it turned out to be inadequate to explain the 
observed high insulating-like resistivity above the transi- 
tion temperature and, even most importantly, the obser- 
vation of CMR in stoichiometric phases like Tl2Mn207»2 
Conversely, the more recent half-metallic ferromagnetic 
model ; 10 ! 11 which is based on the spin-polarized calcula- 
tion of the density of states (DOS) within the Density 
Functional Theory (DFT), provides several clues to the 
underlying processes involved in the CMR phenomena 
and since it does not rely on the mixed-valence (Mn 3+ 
and Mn 4+ ) picture it can explain the observation of CMR 
in stoichiometric phases^ 



Controlling and understanding CMR-HMFs phenom- 
ena in manganites within the DFT framework remains 
a great challenge because of two fundamental obstacles: 
(i) strong exchange-correlation effects and the concur- 
rent orbital/lattice/spin correlations which are not well 
described by conventional DFT methods, and (ii) doping- 
induced structural disorder, which unavoidably limits the 
application of quantum mechanical schemes based on re- 
peated unit cells. The drawbacks of DFT in dealing 
with insulating transition metal oxides can be corrected 
by employing more sophisticated methods such us hy- 
brid functionalsi 2 - which have been proven to provide 
substantially improved structural, electronic and mag- 
netic properties* 1 ^ - — thanks to the inclusion of a portion 
of exact 'non-local' exchange. Beyond-DFT approaches 
are particularly necessary to correctly predict structural 
distortions, magnetic energies and bandgap in LaMn03, 
which are wrongly described by DFT i 13 i 17 i 18 DFT is un- 
able to describe the JT instabilities, and stabilize a metal- 
lic FM solution instead of the experimentally observed 
JT-distorted AFM insulating stated As for the struc- 
tural disorder, the possibility to circumvent the problem 
by realizing CMR-HMF behaviors in stoichiometric sam- 
ples would represent a substantial benefit not only for 
theory, but also for the experimental and technological 
research, thanks to the higher degree of control and ma- 
nipulability of the relevant physical processes. However 
the task is hard, and up-to-date very few stoichiometric 
CMR compounds have been identifiedi^i 1 ^ 

Boosted by the recent experimental observations 
of an insulator-to-metal transition (IMT) in dense 
LaMnC-3^2 at about 32 GPa, in this article we explore 
the structural and magnetoelectric response of LaMnOs 
upon compression up to 150 GPa. 

We first recall the basic properties of LaMnOs. At 
zero pressure and low temperature LaMnOs is a type- 
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TABLE I. Structural parameters of LaMnOs: low-temperature (4.2 K) experimental data^ versus fully optimized PBE and 
HSE results with different values of the exact-exchange mixing parameter a: (PBE), 0.10 (HSE-10), 0.15 (HSE-15), 0.25 
(HSE-25), and 0.35 (HSE-35). Mn-O s , Mn-O m and Mn-O;, represent the short, medium and long Mn-O bondlengths, whereas 
Mn-O m -Mn (°) and Mn-O s / r Mn (°) indicates the corresponding angles. Finally, the JT parameters Q2 and Q3 are defined as: 

Q 2 = 2(1 - s)/y/~(2) and Q 3 = 2(2m - I - s)/y/{6). 
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A AFM insulator (alternating planes of similar spins 
along the c direction) characterized by staggered JT and 
GdFeC>3-type (GFO) distortions, manifested by long (1) 
and short (s) Mn-0 in-plane distances and medium (m) 
Mn-0 vertical ones (JT), and by the tilting of the the 
Mn 3+ 06 octahedra (see Fig. E{d)). These structural in- 
stabilities removes the e g orbital degeneracy and stabilize 
an orthorhombic high-spin (t2g) 3 (e ff ) 1 orbitally ordered 
configuration^ The application of hydrostatic pressure 
progressively quenches the cooperative JT distortions 
and leads to an IMT at P c =32 GPa3 The persistence 
of the structural distortions up to P c indicates that the 
IMT is not a Mott-Hubbard type. This conclusion was 
initially proposed by LDA+U and Dynamical mean field 
theory studie d 17 ' 24 and only very recently was confirmed 
by high pressure Raman measurements ^ Baldini and 
coworkers^ 2 - have also reported the coexistence of do- 
mains of distorted and regular octahedra in the pressure 
range 3-34 GPa, and connected the onset of metallicity 
with the increase of undistorted MnOg octahedra beyond 
a critical threshold. The concomitant presence of two dis- 
tinct phases in this pressure range was confirmed by the 
X-ray absorption spectroscopy experiments of Ramos et 

Our computational study, beside providing a detailed 
microscopic understanding of the IMT and of the as- 
sociated competition between distorted and undistorted 
phases, predicts that the onset of metallicity is associated 
with a FM spin transition and that the FM-metal state 
develops towards a transport half-metal regime at ele- 
vated pressure (« 100 GPa), where the current is nearly 
fully spin polarized. Our results are reported and dis- 
cussed in Scc lIIIl Before that we describe our computa- 
tional setup, which is given in Sec[II] Finally, in Sec lIVI 
we draw a summary. 



II. METHODOLOGY AND COMPUTATIONAL 
ASPECTS 

All calculations were performed using the the Vienna 
Ab initio Simulation Package 2 ^! (VASP) within the 
Heyd, Scuseria, and Ernzcrhof (HSE) hybrid density 
functional scheme^ in which the exchange-correlation 
functional is expressed as a suitable admixture of DFT 
and Hartee-Fock (HF): 

E» S C E = aE^ sr ^+(l-a)E EBE ^+E EBE ^+E BBE 

where /j=0.20A— 1, controls the range separation be- 
tween the short-range (sr) and long-range (lr) part of the 
Coulomb kernel, and a determines the fraction of exact 
HF exchange incorporated. The parameter a, which we 
set equal to 0.15, is chosen so as to provide accurate val- 
ues for the band gap, structural distortions and magnetic 
energies, as discussed below. 

We have used a computational unit cell containing four 
LaMnOa formula units (i.e. 20 atoms) to simulate both 
the Pnma and simple-cubic phase within the ferromag- 
netic (FM) and type-A AFM orderings. Convergence 
tests on the energy cut-off has shown that at low and 
intermediate pressure range the energy difference AE be- 
tween the FM and AFM phases (our most critical quan- 
tity) changes by less than 3 meV per formula unit (f.u.) 
when the energy cut-off is increased from 300 meV to 
400 meV. We have therefore chosen the value 300 meV, 
and made use of a 4x4x4 Monkhorst-Pack k-point grid 
for both the exchange-correlation DFT kernel and the 
portion (a=0.15) of exact HF exchange. 

We have used the following structural optimization 
procedure. For each volume we have fully relaxed the 
cell shape (i.e., lattice parameters a, b, and c as well as 
the corresponding angles between them) and all internal 
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FIG. 1. (Color online) PBE & HSE band structure of 
LaMnOs along the ZrSR (HSE) and LPSRZP (PBE) ob- 
tained using different values of the mixing parameter at the 
corresponding optimized geometry. The Fermi level is set to 
zero. For PBE we have included more high-symmetry lines in 
order to show that both e g bands above and below Ef cross 
the Fermi level. The PBE bands along the ZPSR are shown 
as (red) dashed lines in panel HSE-10. 



structural degrees of freedom (all atomic positions) . This 
complete geometrical optimization allowed us to tread 
the structural path from the P nm a to the cubic phases. 
We have tested about 30 different volumes. The equation 
of states (Pressure- Volume curves) were calculated in two 
different ways: (i) by computing the pressure directly 
from the stress tensor (as automatically done by VASP) 
and (ii) by applying the Birch-Murnaghan equation^ 
Both routes lead to the same result. 



TABLE II. The band gap (A, eV), magnetic moment (m, 
fiB /Mil), and relative energy with respect to the FM ordering 
calculated by HSE for different values of the mixing parameter 



The results, collected in Tab|I] (optimized geometry) and 
Tabim (bandgap, magnetic moment and magnetic ener- 
gies), demonstrate that the best choice is a=0.15. For 
this value of the mixing parameter HSE delivers (see 
TabU and TabHIf (i) an insulating bandgap, A=1.45, 
well within the measured data, (ii) a correct description 
of the critical cooperative Jahn- Teller (JT) parameters 
Q2 and Q3, and (hi) a AFM magnetic ground state in 
agreement with the experimental findings. Smaller mix- 
ing factors (0 and 0.10) lead to a significant underesti- 
mation of the JT parameters Q2 and Q3 (see Tab|TJ and 
to a much too small bandgap (which is actually zero in 
PBE, see Tab HI)) , and to the stabilization of the FM spin 
arrangement (see Tab|n|, in contrast to the experimental 
situation. This is in line with previous conventional DFT 
studies i 18 i 33 ~ — We note that by using the experimental 
structure a small bandgap of about 0.2 eV is opened at 
PBE level, which is still too small as compared to the 
measured level. Conversely, an higher a (0.25 and 0.35) 
correctly favors the AFM ordering but overestimates the 
band gap (2.47 and 3.41 for a=0.25 and a=0.35, respec- 
tively). The value of the gap obtained for cv=0.25 is very 
similar to the corresponding B3LYP value, 2.3 eV^ and 
to the HSE value obtained using the experimental struc- 
ture, 2.25 eV4£ 

In terms of the band dispersions the effect of the mix- 
ing parameter is the progressively larger separation of 
the occupied and unoccupied e g bands below and above 
the Fermi energy with increasing a, as shown in Fig. [TJ 
but qualitatively the overall bonding picture remains un- 
changed. In contrast, in PBE the e g sub-bands cross the 
Fermi Energy and form a spurious metallic solution. 
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A. Choice of the mixing factor a: ground state 
properties of LaMnOs 

In order to determine the HSE mixing factor ev, i.e. the 
fraction of non-local HF exchange included in the hybrid 
exchange-correlation functional, we have performed a set 
of calculations of the ground state structural, electronic 
and magnetic properties of LaMn03 for different values of 
a: (corresponding to a purely PBE setup), 0.10 (HSE- 
10), 0.15 (HSE-15), 0.25 (HSE-25) and 0.35 (HSE-35). 



III. RESULTS AND DISCUSSION 

The progression of the structural properties of com- 
pressed LaMnC>3 computed by HSE as a function of 
v=V/Vq is summarized in Fig. [2J whereas the corre- 
sponding development of the electronic and magnetic 
properties is shown in Figs. |3l [5] and [6] In the pres- 
sure range 0-35 GPa, for which experimental data are 
available, our results are in very good agreement with 
measurements in terms of: (i) the pressure- volume equa- 
tion of states and bulk modulus B (B Bx p*=108 GPa, 
B HSE =104 GPa, see Fig. [^a)), (ii) the pressure- 
induced changes in the structural parameters (Fig. HJb)), 
and (iii) the concurrent suppression of the JT modes 
Q2 and Q3 and the band gap at the same compres- 
sion (i^2=0.82, slightly smaller than the experimental 
one, V/V o =0.86, see Fig. ^b) and Fig. E^a)); the 
P=0 HSE gap opened between occupied and empty e g 
states, E g =1.45 eV (Fig. G3c)), is well within the mea- 
sured range, 1.1-1.7 eVi 31 ' 32 Similarly the HSE ground 
state values of Q2 and Q3 match exactly the experimen- 
tal values^ 

The incremental compression of LaMnOs leads to a 
continuous structural transformation from the P=0 dis- 



4 




FIG. 2. (Color online) Evolution of the structural properties of LaMnOa as a function of u=V/Vo as predicted by HSE 
(triangles) and compared with the experimental data (squares) taken from Refs. I20I . I2H and I30I . (a) Equation of state (the 
inset indicates the value of the bulk modulus); (b) Structural parameters: at ^3=0. 76 LaMnC>3 undergoes an orthorhombic 
(o/i/ \/ c /2) to cubic (a = b = -1/ c/2) transition marked by the vertical line, (c) Progressive quenching of the cooperative 
Jahn- Teller local modes Q2 = 2(1 — s)/y2 and Q3 = 2(2m — / — s)/y/E with increasing pressure, where I, s and m indicate the 
long, short and medium Mn-O bondlength; the JT modes are almost completely quenched at the onset of metallicity, marked 
by the vertical line at W2=0.82 (see Fig. [3j) ; (d) Side (top) and top (bottom) view of the Pnma (left, V/Vo=l) and cubic 
(right, V/Vo<0.76) phases of LaMn03, underlining the suppression of the JT and GFO structural distortion in the perfect 
cubic phase. 



torted Pnma phase to a perfect cubic structure via a 
gradual quenching of the JT modes, the rectification 
of the GFO tilting angles and the alignment of the a, 
b, and c lattice parameters towards the same value, ~ 
5.1 A at w 3 =0.76 as outlined in FigHVd). The e g 
bands around the Fermi energy (Ep) come progressively 
closer until the gap is closed (Figj3jc-f)). Concomi- 
tantly, the unoccupied t2 g bands is pushed down in en- 
ergy and ultimately crosses the Ep at U2=0.82, the on- 
set of metallicity (see Fig. [3jf)). At this critical vol- 
ume HSE predicts a jump in the relative stability be- 
tween the AFM and FM ordering, with the latter be- 
coming the most favorable one by about 90 meV/f.u., 
as illustrated in Fig. [3{b). At low/intermediate com- 
pressions (V/Vo>W2=0.82) the data displayed in Figj3jb) 
shows a strong competition between the AFM and FM 
phases. HSE predicts a crossover between the AFM and 
FM phases at v\=0.Qb (corresponding to a pressure of 
11 GPa), below which the AFM and FM ordering be- 
come almost degenerate (AE<12 meV/f.u.). Consider- 
ing that in the FM phase the JT/GFO distortions are 
almost completely inhibited (see Fig. 0]), this result 
strongly supports the latest Ramao^ and X-ray absorp- 
tion spectroscopy^ studies reporting the formation of a 
mixed state of domains of distorted and regular MnOg 
octahedra in the range 13-34 GPa, which compare well 
with the corresponding theoretical pressure range, 11- 



50 GPa (t>2<V/Vo<u i). The coexistence of distorted 
and undistorted octahedra is clarified in Fig. [4] where 
we plot the comparison between the evolution of the JT 
distortions in the FM and AFM phases (panel (a)), and 
the energy- volume phase diagram (panel (b)). In the 
FM phase the Q2 and Q3 drop down to zero at about 
11 GPa (V/Vo=wi=0.95), indicating that for pressure 
larger than 11 GPa the MnOg octahedra are undistorted. 
Converesly, as already underlined, in the AFM phase the 
octahedra remain distorted until V/Vo=i^2=0.82 (about 
50 GPa). This different behaviour is interpreted graph- 
ically in the insets of Fig. HJa) which represents the 
FM-undistorted (u) and AFM-distortcd (d) octahedra. 
Summing up, at low/high pressures LaMn03 is AFM- 
distorcd/FM- undistorted, whereas in the volume range 
■l>2<V/Vo<wi these to phases coexists in a mixed do- 
mains of distorted and undistorted octahedra (see Fig. 

Bib))- 

The FM transition at V/Vo=0.82 comes right before 
an high spin (HS, S=2) to low spin (LS, S=l) moment 
collapse, which is correlated with the e g and t2 g orbital 
occupations as shown in Fig. [6l under compression the 
Mn 3+ ion retains its P=0 fag)^ {e g Y orbital configura- 
tion down to V/Vo=0.80, with a magnetic moment of 3.7 
(jlb\ further compression yields a rapid reduction of the 
magnetic moment down to 1.7 hb, due the redistribution 
of electrons within the 3d shell which ultimately leads to 
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FIG. 3. (Color online) Evolution of the calculated electronic properties and magnetic ordering of LaMnOa upon pressure, (a) 
Band gap: IMT at W2=0.82 (marked by the vertical full line; the dashed line refers to the corresponding experimental onset). 
Measured gaps are taken from Ref. . 31 ' 32 (b) Energy difference AE between the AFM and FM spin arrangements: AFM/FM 
crossover at ui=0.95 (indicated by a vertical dashed line) and stabilization of the FM state at W2=0.82. (c-1) Changes in the 
e g and ti g density of states around the Fermi level with pressure. The dashed (red) lines refer to the Oxygen p states, whereas 
the thick (red) arrow indicates the transfer of one electron from the e g to the tig sub-bands. 



the low-spin configuration (tig^^^ifigf ■ The HS-to-LS 
collapse starts to develop at about V/Vo=0.8 and is fully 
established at exactly the same volume at which the cu- 
bic phase emerges, t>3=0.76. This HS-orthorhombic to 
LS-cubic transition is also reflected in the DOS (Fig. [3] 
(g-i)), whose evolution from V/V o =0.79 to V/V o =0.73 
clearly indicates the transfer of one electron from the 
e g to the t2 g sub-bands and the subsequent realization 
of a nearly FM half-metallic state with a metallic mi- 
nority t2 g band and a quasi-insulating majority channel 
with a residual density of e g electrons at the bottom of 
the conduction band. In order to clarify further the HS- 
to-LS transition we display in Fig. [5] the partial and 
integrated density of states in a wide energy windows 
associated with the e g and t2 g bands at V/Vo=0.79 and 
V/Vo=0.73, which show the transfer of one electron from 
the e g spin up channel at V/Vo=0.79 to the t^g spin-down 
channel at V/V =0.73, which yields to the (t 29 ) ttt (e 9 ) t 
to (t2g)^^{e g )° transition. 

Following the classical work of Nadgorny et al*& we 
have analyzed the spin polarization P n associated to this 
high pressure nearly HM-FM state in order to acquire 
information on the spin-dependent transport properties, 
using the formula suggested by Mazing 

n N t {E F )v Ft + N^Ef^ [ ' 

where N^(Ep), N±(Ef) and VFf, v F , represent the ma- 



jority and minority spin DOS and Fermi velocities, re- 
spectively, and the index n refers to the different types 
of spin polarizations detected in spin-resolved photoc- 
mission measurements (n = 0), and in ballistic (n = 1) 
and diffusive (n = 2) transport experiments. We have 
computed the Fermi velocities by interfacing the VASP 
with the BoltzTrap code,^ and obtained for V/V o =0.70: 
P = 87%, Pi = 80.5%, and P 2 = 71%, and for 
V/V o =0.65: P = 92%, P 1 = 87% and P 2 = 80%, values 
very similar to those reported for the doped CMR man- 
ganite Lao.7Sro.3Mn03.— We can thus conclude, that the 
high pressure FM cubic phase of LaMnOa is a transport 
half-metal. For denser phases (V/Vo<0.65) the majority 
spin band gap (from the lower laying filled t2 g and the 
unoccupied e g band) is progressively reduced to zero at 
V/V o <0.53 (P>300 GPa, see FigE^a)). Being the FM- 
tHM regime the crucial common ingredient of all CMR 
manganites, its realization in the undoped (stoichiomet- 
ric) phase of the CMR parent compound LaMnOs in a 
wide interval of compressions, could help in achieving 
new fundamental insights into the elusive phenomena of 
CMR. 



IV. SUMMARY 

Summarizing, we have traced the development of 
LaMnOs upon pressure and determined a sequence of 
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FIG. 4. (Color online) (a) Development of the JT quantities Qi and Q2 upon pressure for the FM and AFM ordering, (b) 
Energy- Volume curve showing the coexistence of JT-distorted and undistorted phases in the volume range W2iV/Vojwi. This 
volume region is delimited by the vertical lines. In both panels the insets represent FM-undistorted (u) and AFM-distorted (d) 
MnOe octahedra. 
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FIG. 5. (Color online) Partial and integrated density of states 
of LaMn0 3 at V/Vo=0.79 and V/V =0.73 representing the 
pressure induced High-Spin to Low-Spin transition and as- 
sociated e g t (V/Vo=0.79)-^t2 S ; (V/V =0.73) electron trans- 
fer. The arrows indicate the HS and LS configuration at 
V/Vo=0.79 and V/Vo=0.73, respectively, whereas the small 
numbers represent the number of electrons for the specific Im- 
channel (those in squares are those involved in the electron 
transfer). 



highly interconnected structural, electronic, and mag- 
netic phase transitions: (i) At ambient conditions 
LaMnC>3 posses a distorted AFM insulating state. At 
??i=0.95 a competition between (distorted) AFM and 
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FIG. 6. (Color online) (a) Pressure driven magnetic moment 
collapse and corresponding correlation with (b) the e 9 and 
t2 g orbital occupancy. The insets in panel (a) illustrate the 
schematic diagrams of the HS and LS states at low (right) and 
high (left) pressure. The square indicates the experimental 
magnetic moment at ambient pressure taken from Ref.i^ A 
similar transition has been observed in dense MnO.— 



(undistorted) FM configuration begin to evolve, (ii) At 
the critical threshold ^2=0. 82 LaMnC>3 undergoes a non- 
Mott IMT associated with a significant reduction of the 
JT/GFO instabilities, (iii) For i; 3 =0.76 all residual lat- 
tice distortions are suppressed and a perfect cubic phase 
emerges, (iv) At ^4=0.70 the strong crystal field splitting 
between ti g and e g drives a magnetic moment collapse 
from an HS to LS tHM-FM state, manifested by a signif- 
icant disproportionation in the spin-dependent N(Ep), 
Vp, and consequentially P n . Our results, thus, predict 
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that it is possible to establish a tHM-FM regime in a 
stoichiometric manganite at pressure accessible by high- 
pressure technology. 

ACKNOWLEDGEMENTS 

Support by European Community (FP7 grant 
ATHENA), the Chinese Academy of Science (CAS Fel- 



lowship for Young International Scientists) and the 
National Science Foundation of China (NSFC Grand 
51050110444) is gratefully acknowledged. X.-Q.C is 
grateful for supports from the CAS Hundred Talents 
Project and from NSFC of China (Grand Numbers: 
51074151 and 51174188). All calculations have been per- 
formed on the Vienna Scientific Cluster. 



* Corresponding author: cesare.franchmi@univie.ac. at 

1 R. A. de Groot, F. M. Mueller, P. G. van Engen, and K. 
H. J. Buschow, Phys. Rev. Lett. 50 2024 (1983). 

2 M. I. Katsnelson, V. Yu. Irkhin, L. Chioncel, A. I. Licht- 
enstein, R. A. de Groot, Rev. Mod. Phys. 80, 315 (2008). 

3 I. Zutic , J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

4 G. A. Prinz, Science 282 1660 (1998). 

5 J.M.D. Coey, M. Viret and S. von Molnar, Adv. in Phys. 
48, 167 (1999). 

6 G. H. Jonker and J. H. van Santen, Physica 16, 337 (1950). 

7 M. B. Salamon and M. Jaime, Rev. Mod. Phys. 73, 583 
(2001). 

8 C. Zener, Phys. Rev. 82, 403 (1951). 

9 Y. Shimakawa, Y. Kubo, and T. Manako, Nature 379, 53 
(1996). 

10 B. Nadgorny, I. I. Mazin, M. Osofsky, R. J. Soulen, Jr., 
P. Broussard, R. M. Stroud, D. J. Singh, V. G. Harris, 
A. Arsenov, and Ya. Mukovskii, Phys. Rev. B 63, 184433 
(2001). 

11 W. E. Pickett and D. J. Singh, Phys. Rev. B 53, 1146 
(1996). 

12 A.D. Becke, J. Chem. Phys. 98 1372 (1993). 

13 D. Munoz, N. M. Harrison, and F. Illas, Phys. Rev. B 69, 
085115 (2004). 

14 T. Archer, CD. Pemmaraju, S. Sanvito, C Franchini, J. 
He, A. Filippetti, P. Delugas, D. Puggioni, V. Fiorentini, 
R. Tiwari and P. Majumdar, Phys. Rev. B 84, 115114 
(2011). 

15 C Franchini, A. Sanna, M. Marsman, G. Kresse, Phys. 
Rev. B 81, 085213 (2010). 

16 C Franchini, T. Archer, Jiangang He, Xing-Qhi Chen, A. 
Filippetti, S. Sanvito, Phys. Rev. B 83, 220402(R) (2011). 

17 G. Trimarchi and N. Binggeli, Phys. Rev. B 71, 035101 
(2005). 

18 T. Hashimoto, S. Ishibashi, ann K. Terakura, Phys. Rev. 
B 82, 045124 (2010). 

19 S. Ishiwata, T. Nakano, I. Terasaki, H. Nakao, Y. Mu- 
rakami, Y. Uwatoko, and M. Takano, Phys. Rev. B 83, 
020401(R) (2011). 

20 I. Loa, P. Adler, A. Grzechnik, K. Syassen, U. Schwarz, 
M. Hanfland, G. Kh. Rozenberg, P. Gorodetsky, and M. P. 
Pasternak, Phys. Rev. Lett. 87 125501 (2001). 



21 L. Pinsard-Gaudart, J. Rodrii'guez-Carvajal, A. Daoud- 
Aladine, I. Goncharenko, M. Medarde, R. I. Smith, and 

A. Revcolevschi, Phys. Rev. B 64, 064426 (2001). 

22 M. Baldini, V. V. Struzhkin, A. F. Goncharov, P. Pos- 
torino, and W. L. Mao, Phys. Rev. Lett. 106 066402 
(2011). 

23 Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, I. Koyama, 
M. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
and Y. Endoh, Phys. Rev. Lett. 81, 582 (1998). 

24 A. Yamasaki, M. Feldbacher, Y.F. Yang, O.K. Andersen, 
K. Held, Phys. Rev. Lett. 96, 166401 (2006). 

25 AY. Ramos, N. M. Souza-Neto, H. C N. Tolentino, O. 
Bunau, Y. Joly, S. Grenier, J.-P. Itie, A.-M. Flank, P. La- 
garde and A. Caneiro, Europhys. Lett. 96 36002 (2011). 

26 G. Kresse and J. Furthmiiller, Comput. Mat. Sci. 6, 15 
(1996). 

27 G. Kresse and D. Joubert, Phys. Rev. B, 59, 1758 (1999). 

28 J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 
118, 8207 (2003); Erratum: J. Chem. Phys. 124, 219906 
(2006). 

29 Francis Birch, Phys. Rev. 71 809 (1947). 

30 J.B.A.A. Elemans, B. van Laar, K.R. van der Veen, and 

B. O. Loopstra, J. Phys. Chen. Solids 3, 238 (1971). 

31 T. Arima, Y. Tokura, and J.B. Torrance, Phys. Rev. B 48 
17006 (1993). 

32 T. Saitoh, A. E. Bocquet, T. Mizokawa, H. Namatame, A. 
Fujimori, M. Abbate, Y. Takura, M. Takano, Phys. Rev. 
B 51, 13942 (1995). 

33 H. Sawada, Y. Morikawa and K. Terakura, N. Hamada, 
Phys. Rev. B 56, 12154 (1997). 

34 H. Sawada and K. Terakura, Phys. Rev. B 58, 6831 (1998). 

35 E. A. Kotomin, R. A. Evarestov, Yu. A. Mastrikov and J. 
Maier, Phys. Chem. Chem. Phys. 7, 2346 (2005). 

36 C Franchini, R. Kovacik, M. Marsman, S. Sathyanarayana 
Murthy, J. He, C Ederer, and G. Kresse. larXlv:llll. 15281 
(2011). 

37 J. Kunes, Alexey V. Lukoyanov, Vladimir I. Anisimov, 
Richard T. Scalettar and Warren E. Pickett, Nat. Mat. 
7, 198 (2008). 

38 I.I. Mazin, Phys. Rev. Lett. 83, 1427 (1999). 

39 G.K.H. Madsen and D.J. Singh, Comput. Phys. Commun. 
67 175 (2006). 



